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Abstract — Background: The accumulation of deleterious mu- 
tations of a population directly contributes to the fate as to 
how long the population would exist. Muller's ratchet provides 
a quantitative framework to study the effect of accumulation. 
Adaptive landscape as a powerful concept in system biology 
provides a handle to describe complex and rare biological events. 
In this article we study the evolutionary process of a population 
exposed to Muller's ratchet from the new viewpoint of adaptive 
landscape which allows us estimate the single click of the ratchet 
starting with an intuitive understanding. 

Methods: We describe how Wright-Fisher process maps to 
Muller's ratchet. We analytically construct adaptive landscape 
from general diffusion equation. It shows that the construction 
is dynamical and the adaptive landscape is independent of the 
existence and normalization of the stationary distribution. We 
generalize the application of diffusion model from adaptive 
landscape viewpoint. 

Results: We develop a novel method to describe the dynamical 
behavior of the population exposed to Muller's ratchet, and 
analytically derive the decaying time of the fittest class of 
populations as a mean first passage time. Most importantly, 
we describe the absorption phenomenon by adaptive landscape, 
where the stationary distribution is non-normalizable. These 
results suggest the method may be used to understand the 
mechanism of populations evolution and describe the biological 
processes quantitatively. 

Index Terms — Wright-Fisher process, adaptive landscape, sta- 
tionary distribution, mean first passage time 

Background 

Muller's ratchet is a mechanism that has been suggested 
as an explanation for the evolution of sex (T). For asexually 
reproducing population, without recombination, chromosomes 
are directly passed down to offsprings, as a consequence 
the deleterious mutations accumulate so that the fittest class 
loses. For sexually reproducing population, because of the 
existence of recombination between parental genomes, a 
parent carrying high mutational loads can have offspring 
with fewer deleterious mutations. The high cost of sexual 
reproduction is thus offset by the benefits of inhibiting the 
ratchet Q. Muller predicted ratchet mechanism in 1964. 
Muller's ratchet has received growing attention recently. Most 
studies of Muller's ratchet are related to two issues. First 
is that without recombination, the genetic uniformity of the 
offspring leads to much lower genetic diversity, which is 



likely to make it more difficult to adapt to [3]. Second is that 
population lacking genetic repair should decay with time, due 
to successive loss of the fittest individuals HO. In addition, 
Muller's ratchet is relevant to some replicators, endosymbionts 
and mitochondria. In order to assess the relevance of Muller's 
ratchet, it is necessary to determine the rate (or the time) for 
the accumulation of deleterious mutations [6|. It is widely 
recognized that the rate of deleterious mutations being much 
higher than that of either reverse or beneficial mutations 
can be a serious threat to the survival of populations at the 
molecular level J4). It would have to be taken into account in 
any quantitative theory of the ratchet mechanism Q. 

Some biologists also have noticed a framework was 
needed to construct such as [7|[6|. The potential evolutionary 
importance of Muller's ratchet mechanism would make it 
desirable to carry out careful quantitative studies on its 
operation 0. There has been a continuous progress in 
theory. The simplest and earliest mathematical model is 
the pioneering work in (8). The author In J5] described 
the evolutionary process of asexual population of constant 
size on the condition of deterministic mutation selection 
balance according to Wright-Fisher dynamics. And the 
author found numerical evidence of relation between the 
total number of individuals and the average time between 
clicks of the ratchet. Authors in |9| treated pioneering model 
as diffusion approximation, they produced more accurate 
predictions over the relatively slow regime. They noted the 
increasing importance of selection coefficient for the rate 
of the ratchet for increasing values of the total number of 
individuals. The author in iflOl employed simulation and 
analytical approaches to Muller's ratchet and estimated 
how different between the distribution of mutations within 
a population and a Poisson distribution. The author in 
121 obtained diffusion approximations for three different 
parameter regimes, depending on the speed of the ratchet. 
The model mathematically embodied as stochastic differential 
equation and shed new light on (9). Authors in also 
mapped Muller's ratchet to Wright-Fisher process, and got 
the prediction of the rate of accumulation of deleterious 
mutations when parameters lie in the fast and slow regimes 
of operation of the ratchet. In the present article, inspired by 



recent works such as (6), we treat the evolutionary process as 
Wright-Fisher model for diffusion approximation, analytically 
construct adaptive landscape, where the stationary distribution 
is non-normalizable. This method expands the application 
of the model in [6|. The concept of adaptive landscape 
is proposed by Sewall Wright to quantitatively describe 
the complex biological phenomena as ifTTl . Here adaptive 
landscape is clearly quantified as potential function in [12| 
from the physical point of view. 

The key concept in constructing adaptive landscape is of 
potential function as a scalar function. There is a long history 
of definition, interpretation, and generalization of potential. 
Potential has also been applied to biological systems in 
various ways. The usefulness of a potential reemerges in the 
current study of dynamics of gene regulatory networks Q~3], 
such as its application in genetic switch |14||15]|16]. The 
role of potential is the same as that of adaptive landscape. In 
this article, we do not distinct them. 

In this article we present an approach to Muller's ratchet, 
mapping a single click of Muller's ratchet to one dimensional 
Wright-Fisher model where two allele classes are used to 
depict all different types of individuals in a population. 
The Wright-Fisher model predicts the rate of Muller's 
ratchet and suggests the accurate rate of accumulation of 
deleterious mutations could be derived more quickly by 
adaptive landscape. In addition we will demonstrate this 
method is independent of the existence and normalization of 
stationary distribution. Let we further address one dimensional 
Wright-Fisher process to make our points more clearer. 

Methods 

Model description 

We consider here an important and widely applied process- 
Muller's ratchet which is mapped to Wright-Fisher process. 
Muller's ratchet is the process by which genomes of a finite 
population composed of asexual individuals accumulate dele- 
terious mutations in an irreversible manner [17||18|. Consider 
a population of haploid asexual individuals with discrete 

generations t = 0, 1, 2, The starting point in a generation 

is regarded as adult stage, after all selection has occurred and 
immediately prior to reproduction. New mutations occur at 
reproduction and all mutations are assumed to deleteriously 
affect viability but have no effect on fertility. Supposed 
population size is always fixed each generation, then there 
are always alleles with constant number in the gene pool. 
Here we consider one locus with two alleles A and a, that 
is, there are two classes in the haploid asexual population. 
One class is individuals with allele A while the other is 
individuals with allele a. Supposed mutation from allele A 
to a is deleterious, we assume population of constant size 
N > 1, generations are non-overlapping. The lifecycle of the 
individuals in the population is from adults to juveniles, during 
which we consider the change of allele A in the presence of 
irreversible mutations, selection and random drift. The change 
of fittest allele frequency is exposed to Muller's ratchet. Given 



above assumption, we can study the evolution of allele A. The 
frequency of allele A in generation t is x t while that of allele 
a is 1 — x t , x t = 0, 1/N, . . . , 1. Let p be the probability that 
an offspring of an adult with allele A is an individual with 
allele a. It is labeled by Mi o, that is, M±fi — p. Analogously, 
Mo .o = 1 — M> Mo i = 0, Mm = 1. The relative viability of 
individuals with allele A is v = 1 while that of individuals 
with allele a is v\ = 1 — er, where a can be treated as effective 
selection rates associated with deleterious mutations. Then in 
generation (t + 1), the frequency of allele A is 

(1 - p)x t 



X t +1 = 



(1) 



1 — a + ct(1 — p)x t 

Under the general diffusion approximation, frequency Xt is 
treated as continuous quantities x, and this also leads to the 
distribution of the frequency of allele A being probability 
density. Let p = p(x, t) be the probability density of the 
frequency of allele A being x at time t. The diffusion equation 
obeys the following formula 1 19l [|20ll 



d t p = dl[V{x)p/2]-d x [M{x)p], 



with 



M(x) = 
V{x) = 



x[(o — fi) — a(l — p)x] 

1 — a + cr(l — p)x 
x(l — x) 
N ' 



(2) 



(3) 



(4) 



where M(x) represents the average transition rate of x or 
driving force [19] [21] and V(x) is the corresponding variance. 

Adaptive landscape 

We can also depict the same evolutionary process by the 
symmetric equation 

d t p = d x [eD(x)d x -f(x)]p (5) 

with 

x[(a — p) — cr(l — p)x] 1 — 2x 

1 — a + cr(l — p)x 
x(l — x) 



2N 



eD{x) 



2N 



(6) 
(7) 



(8) 



Adaptive landscape is directly given by 

We are interested in the dynamical property of adaptive 
landscape, so we treat $ and <£>/e no different in this respect, 
that is, for convenience we can take e = 1 of eD(x). 



fix) 
eD(x) 
2Np(l - a) 
1 — ap 
2N(l - p) 



= lMt) dX (9) 

2/V//Y1 - rr^ 

ln(l — x) — lnx(l — x) 

ln(l — a + x<j{1 — p)). 
1 — a p 

Here the adaptive landscape is composed of three terms. The 
first term and the third term quantify the effect of irreversible 



mutation and selection respectively, the second term quantifies 
the effect of random drift. 

In addition, the symmetric Eq.(5) has two advantages. On 
the one hand, the adaptive landscape is directly read out 
when the detailed balance is satisfied. On the other hand, the 
constructive method is dynamical, independent of existence 
and normalization of stationary distribution. We call f(x) 
directional transition rate, integrating the effects of M{x) and 
the derivative of V(x). Directional transition rate can give 
equilibrium states when f(x) is linear form. 

Under the very natural boundary condition satisfying the 
probability flux of the system is zero, we assume the probabil- 
ity flows in x £ [0,1]. The stationary distribution for diffusion 
approximation is given by 



p{x) 



1 



: exp 



$(a?) 



The stationary distribution can also be expressed as 

' ^-i-^Qx-i) 

1 — <T/i 

2N(l - u) , , , „\ 

— ^ ^ln (1-a + xa (1-p)) , 

1 - (T/i J 



- lnx - 



(10) 



it has the form of Boltzmman-Gibbs distribution 11221 . so the 
scalar function $(a;) naturally acquires the meaning of poten- 
tial energy lfl3l . The value of Z determines the normalization 
of p from the perspective of probability, and the finite value of 
Z manifests the normalization of p. The constant e holds the 
same position as temperature of Boltzmman-Gibbs distribution 
in statistical mechanics. But it does not hold the nature of 
temperature in Boltzmman-Gibbs distribution. 

Results and Discussion 

Previous works mainly focus on the parameters ranges that 
either lie at the lower mutation rates regime or the higher mu- 
tation rates regime based on stochastic differential equations. 
We analyze the evolutionary process in all parameters regimes 
based on Fokker-Planck equation by adaptive landscape. We 
investigate the single click time analytically. 

Irreversible mutation, selection and random drift balance 

To understand the mechanism of Muller's ratchet, a full 
characterization of dynamical process is a prerequisite for 
obtaining more accurate decaying time. Here we study it by 
adaptive landscape in detail. 

Concretely we divide the regime of mutation rates into 
three sections. Mutation rates p £ (0, 1/(2N — 1)], selection 
rates a £ (p,2p/(l + p)}, this condition corresponds to 
the case that selection rates and mutation rates lie in the 
lower regime. In the medium regime mutation rates p £ 
(1/(2N - 1), l/(2N - 2)], selection rates a £ (p, 1). In the 
regime of high mutation rates, it is possible that the number 
of mutation is greater than one in the process. Concretely, 
mutation rates p £ (l/(2iV — 2), 1), selection rates a £ (p, 1). 
The adaptive landscape of full parameter regimes is visualized 



as Fig.l, From the expression and visualization of adaptive 
landscape &(x), we may find there are two singular points 
and 1 of adaptive landscape, characterized by infinity value 
in Fig.l. In the low and medium mutation rates regimes, 
singularity indicates the population stable, but in the high 
mutation rates regime, singularity from x = indicates the 
population stable while singularity from x = 1 means the 
population unstable. Fig.l demonstrates the whole process of 
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Fig. 1. Adaptive landscape against allele frequency x with mutation 
rates and selection rates in all regimes, x label represents allele 
frequency, p label represents mutation rates while vertical label is 
corresponding value of $. Two stable states(maximums of adap- 
tive landscape) occur in the low and high mutation rates regimes. 
Assume population size is constant 7V=50. In the lower regime, 
green curve represents £t=0.000005, (7=0.00005, blue curve stands 
for /x=0.000005, cr=0.019, red curve represents ^=0.01,rr=0.019. In 
the medium mutation rates regime, green curve represents /i=0.01015, 
<7=0.015, blue curve stands for /i=0.01015, <r=0.05, red curve inves- 
tigates /i=0.01015, o=0.9. In the higher mutation rates regime, green 
curve represents p=0.4, a=0.9, blue curve stands for p=0.02, rr=0.1, 
red curve describes the case for p=0.02, <r=0.05. 

the population evolution including the forming and losing the 
fittest class. With increasing selection rates the fittest class A 
is forming quickly while with increasing mutation rates the 
fittest class A loses easily. In the lower mutation rates regime 
it describes the population is likely to move to the fittest class 
with increasing selection rates, the process is dominated by 
selection. Blue and red curves manifest the losing process of 
allele A with accumulation of irreversible mutations. Because 
the mutation rates are lower, selection rates dependent of 
mutation rates are lower, these factors result in the difficult 
change of fittest class. As a consequence there are two fittest 
classes in the process. In the middle mutation rates regime, it 
demonstrates the transition process between lower and higher 
mutation rates. The number of mutations is also fewer than 
one in one generation in the process. In the higher mutation 
rates regime, green curve describes the population is likely 
to move to the fittest class so that the population exists in the 
form of coexistence of A and a. Red and blue curves manifest 
the losing process with accumulation of irreversible mutations. 
Because mutation rates are higher, as a consequence allele a 



occurs to a certian percentage. Selection rates dependent of 
mutation rates tend to survive allele A. There are two stable 
states in the process under the balance of irreversible mutation 
and selection. The evolutionary process is dominated by the 
irreversible mutations, the fittest class A loses. So we can draw 
the conclusion that the click process occurs when there are two 
stable states in the process. 

The single click time for Muller's ratchet 

We visualize adaptive landscape, then one may wonder 
about how the population moves from one peak to another 
and how long it might be to move from one maximum 
to another. The process was first visualized by Wright in 
1932. In addition, the problem of transition from metastable 
states is ubiquitous in almost all scientific areas. Most of 
previous works encounter finite potential barriers from the 
physical point of view. Interesting issue here is that we touch 
upon infinite potential barriers under the circumstance of well 
defined two stable states. Then we manifest the derivation 
of a single click time. The time of a click of the ratchet is 
recognized as the random time of loss of the fittest class as 
0. 

The single click time is well defined when two fittest classes 
exist in the process. It means the interval time from the 
forming of one fittest class to its extinction. The corresponding 
processes are that there are two well-defined stable states in 
the system. Corresponding graphs of adaptive landscape is the 
following: 

To evaluate the single click time and show the further power 
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Fig. 2. Adaptive landscape with two stable states in the low and high 
mutation rates. 7V=50, blue curve represents /i=0.02, (7=0.1 while 
green curve stands for ^ = 0.000005, <r=0.00005. 

of adaptive landscape, in the following we will demonstrate 
how the single click time from one stable state to another is 
derived in this framework. After straightforward calculation, 
backward Fokker- Planck corresponding to Eq.(5) can be 
expressed with the property of time homogeneous in the 
following form E31E41 

d t p = (f{x) + eD'(x))d x p + eD(x)dlp, (11) 



where the prime means the derivative of D(x). General single 
click time T(x) dependent on initial Dirac function satisfies 



if(x) + eD'{x))d x T{x) + eD{x)d 2 x T(x) 
With boundary condition 

T'(l) = 0,T(0) = 0. 



(12) 



(13) 



Because of irreversible mutations from A to a, the boundary 
is reflecting at x = 1, the boundary is absorbing at x = in 
the absence of mutation. 

The general solution corresponding to Eq.(12) is 



T(x) 



1 



eD(y) e M-®(y))dy 
J x f(x')/eD(x')dx'(e 



/y 
exp($(z))riz, (14) 



here $(x) = / f(x')/eD(x')dx'(e = 1). 

Here the evolutionary process naturally occurs when x g 
[0, 1]. In the process, there are two important states x* , x^. 
Interval (0,1) contains an potential well at x* , f(x*) = 0, 
f'(x*) < 0, and an potential barrier at x* s , f{x* s ) = 0, 
f'(x^) > 0. Then one can use asymptotic estimation to 
simplify Eq.(14) as 

Ti-s-o = lim T(x) 

x— >1 

27rexp($(a;5) - ®{x*)) 
D(x*)^<5>"(x*)<&"(x%) 

dx. 



/(*) 



(15) 



The second term is for downhill relaxation. The first term is 
for barrier crossing. This method for mean first passage time 
is consistent with l25l . 

In the high mutation rates regime, where Xg approximates to 
a stable state which is near enough to 1, x* corresponds to the 
saddle point that the population lies between the two stable 
states and x$. The single click time approximates to 



Ti^q = lim T(x) 

X— ¥\ 



= lim 



1 



, x exp(-$(y))dy / exp($(z))dz 

x^ij eD(y) J y 



lim 2N 

x-yl 



2JVm(ct-1) 

(l-y) 



(1 — <7 + ycr(l - fi)) i-"" 

L (l - Z) ^ 



2W(l-n) 

x(l — a + az(l — pi)) 1 -"i t dz 
27rexp($(a;5) - $(x*)) 
D{x*)y/<S>"{x*)$"{xl) 



N 



,5-2N 



(4- (7 + 3(7^) 



O(l) 



(16) 



3 (2 - cr/x + 3ct) 

Here x'^ and x* are the zero points of f(x) respectively. For 
the lower mutation rates regime, the single click time can be 
calculated similarly, where X3 approximates to 1, it can be 



derived by mean field method, x* corresponds to the saddle 
point that the population lies at the lowest potential, x* is the 
zero point of f(x). 

Analogous to the derivation of T^q, we can calculate 



0-M 



= lim 

K->0 



eD(y) 



lim 2N 





exp(-$(y))dy / exp($(z))dz 



2JV M (<7- 1) 

(1-y) 



1 (1 - a + ya(l - fi)) 



/ z'^l-z)' 
>y 



7fi-l — 2Nii(cr-l) 



x(l — a + az(l — (J,)) dz 



(17) 



And To-j.1 goes to infinity as a polynomial way. Here the 
single click time asymptotically conforms to exponential form 
according to Eq.(15), though the shapes of adaptive landscape 
are different. And once the population arrives at the stable 
state x = 0, the population is absorbed as Eq.(17) describes. 
The accurate graph of single click time Ti^o is the following: 
Fig. 3 describes the single click time increases with increasing 
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Fig. 3. Exact numerical results about the single click time corre- 
sponding to Eq.(15). Horizontal coordinates represent selection rates 
a and mutation rates /i, vertical coordinate corresponds to the value 
of the click time Ti^o- 

mutation rates and decreases with increasing selection rates. 
This is consistent with adaptive landscape in the lower and 
higher regimes in Fig.l. Asymptotic results corresponds to 
formula (16). 

The estimated results are more accurate when the parame- 
ters lie in the high mutation rates regime. 

Discussion 

The present article presents a different approach to es- 
timate the single click time of Muller's ratchet. Previous 
works mainly based on the stochastic differential equations 
for diffusion approximation to study Muller's ratchet such as 
ifTUl El. This method can result in the complex of calculating 



stochastic differential equations. Inspired by a recent work as 
ll6l , we connect it to one locus Wright-Fisher model with 
asexual population including N haploid individuals. Direct 
classical Wright-Fisher model for diffusion approximation 
can reduce complex calculation to solve matrix equations, 
this method is functional especially when the dimension of 
the matrix is higher. Our theoretical results generalize the 
application of the model proposed in |6|. And it does not 
need the existence and normalization of stationary distribution. 
Our method investigates the global dynamical property of the 
system more directly. In addition, these results demonstrate the 
single click time is approximately exponentially distributed. 
Most importantly, we depict the absorption phenomena from 
the adaptive landscape viewpoint. 

To summarize, we have obtained two main sets of results in 
the present work. First, we demonstrate a transformation of a 
simple Fokker-Planck equation so that the adaptive landscape 
can be explicitly read out as a potential function. Such demon- 
stration suggests this concept applicable to all levels of systems 
biology. This not only allows computing click time of Muller's 
ratchet straightforward, but also makes the stability analysis 
of the ratchet system intuitive. Second, in this framework the 
derivation of single click time is simple even for situations that 
steady state distribution is nonnormalizable. Hence it allows us 
another way to handle absorbing boundary condition. In this 
perspective our work may be a starting point for estimating 
the click time for Muller's ratchet in more general situations. 

Acknowledgments 

We would like to thank Song Xu, Pengyao Jiang, Fangshu 
Cui, Quan Liu, Wei Zheng and others in our lab for many 
conversations. This work was supported in part by the National 
973 Projects No. 2007CB914700 and No. 2010CB529200 
(P.A.)and in part by the project of Xinyang Normal university 
No. 20100073(S.J.). 

References 

[1] J. Maynard Smith. The Evolution of Sex. Cambridge University Press, 
1978. 

[2] A. Etheridge, P. Pfaffelhuber, and A. Wakolbinger. How often does 
the ratchet click? Facts, heuristics, asymptotics. Trends in stochastic 
analysis, pages 365-390, 2009. 

[3] K.P. Lampert and M. Schartl. A little bit is better than nothing: 
the incomplete parthenogenesis of salamanders, frogs and fish. BMC 
biology, 8(1):78, 2010. 

[4] L.P Maia. Analytical results on Muller ratchet effect in growing 
populations. Physical Review E, 79(3):32903, 2009. 

[5] N.H. Barton. Mutation and the evolution of recombination. Philo- 
sophical Transactions of the Royal Society B: Biological Science, 
365(1552):1281-1294, 2010. 

[6] D. Waxman and L. Loewe. A stochastic model for a single click of 
Muller's ratchet. Journal of theoretical biology, 264(4): 1120-1 132, 
2010. 

[7] J. Felsenstein. The evolutionary advantage of recombination. Genetics, 

78(2):737-756, 1974. 
[8] J. Haigh. The accumulation of deleterious genes in a population-muller's 

ratchet. Theoretical Population Biology, 14(2):25 1-267, 1978. 
[9] W. Stephan, L. Chao, and J.G. Smale. The advance of Muller's ratchet in 

a haploid asexual population: approximate solutions based on diffusion 

theory. Genetics Research, 61(03):225-231, 1993. 
[10] D.D.G. Gessler. The constraints of finite size in asexual populations and 

the rate of the ratchet. Genetics Research, 66(03):241-253, 1995. 



[11] H.P. de Vladar and N.H. Barton. The contribution of statistical physics 
to evolutionary biology. Trends in Ecology and Evolution, 26:424—432, 
2011. 

[12] P. Ao. Laws in darwinian evolutionary theory. Physics of Life Review, 
2:117-156,2005. 

[13] P. Ao. Potential in stochatic differential equation: novel construction. 

Journal of physics A: mathematical and general, 37:25-30, 2004. 
[14] X. M. Zhu, L.Yin, L. Hood, and P. Ao. Robustness, stability and 

efficiency of phage lambda genetic switch: dynamical structure analysis. 

Journal of bioinformatics and computational biology, 2(4):785-817, 

2004. 

[15] X. M. Zhu, L.Yin, L.Hood, and P. Ao. Calculating biological behaviors 
of epigenetic states in the phage lambda life cycle. Functional & 
integrative genomics, 4(3): 1 88-195, 2004. 

[16] J. Liang and H. Qian. Computational cellular dynamics based on the 
chemical master equation: A challange for understanding complexity. 
Journal of Computer Science and Technology, 25:154-168, 2010. 

[17] HJ. Muller. Some genetic aspects of sex. American Naturalist, 
66(703):118-138, 1932. 

[18] H.J. Muller. The relation of recombination to mutational advance. Muta- 
tion Research/Fundamental and Molecular Mechanisms of Mutagenesis, 
l(l):2-9, 1964. 

[19] M. Kimura. Diffusion models in populayion genetics. Journal of Applied 

Probability, 1:177-232, 1964. 
[20] WJ. Ewens. Mathematical Population Genetics. Cornell University, 

2004. 

[21] J.H. Gillespie. Population Genetics: a Concise Guide. The Johns 
Hopkins University press, 1998. 

[22] P. Ao. Emerging of stochastic dynamical equalities and steady state ther- 
modynamics from Darwinian dynamics. Communications in Theoretical 
Physics, 49(5):1073-1090, 2008. 

[23] van Kampen N.G. Stochatic Processes in Physics and Chemistry. 
Amsterdam, 1992. 

[24] B. 0ksendal. Stochatic Differential Equations: an Introduction with 

Applications. Springer, 2003. 
[25] H. Qian. Nonlinear stochastic dynamics of mesoscopic homogeneous 

biochemical reaction systemslan analytical theory. Nonlinearity, 24:19- 

49,2011. 



